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ABSTRACT 

This document is intended to acquaint users of the Goddard Trajec- 
tory Determination System (GTDS) Solar/Lunar/Planetary (SLP) 
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with Appendices B through D, constitute a user’s manual for the 
GTDS SLP Ephemeris Files. 
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SECTION 1 


INTRODUCTION 


Using a Jet Propulsion Laboratory (JPL) planetary ephemeris tape as a data 
source, the Goddard Trajectory Determination System (GTDS) has the capability 
of generating Solar /Lunar/Planetary (SLP) ephemeris files which contain data 
in the form of Chebyshev polynomial coefficients. 

The current JPL Export Ephemeris (DE-19) contains the most accurate predic- 
tions (in terms of a long-span time interval) of lunar and planetary motion avail- 
able for the time interval from December 30, 1949 to January 5, 2000. It con- 
sists of three tapes which collectively span the interval as follows: 

• December 30, 1949 to December 29, 1969 

• November 19, 1969 to February 22, 1984 

• January 13, 1984 to January 5, 2000. 

Depending upon the particular body represented, the stepsize for the ephemeris 
data contained on these tapes is either 1/2 day, 2 days, or 4 days. This tabular 
ephemeris data can be used directly via use of the JPL software which provides 
interpolated values of position and velocity vectors of any requested set of bodies 
relative to any requested central body. (A more detailed description of the JPL 
Ephemeris Tape System is given in Appendix A. ) 

However, more efficient use of the tabular data is provided by the GTDS file 
generation procedure by utilizing Chebyshev polynomial curve-fitting techniques 
and by allowing the SLP ephemeris representation of concern to be given in terms 
of long arc lengths in the neighborhood of a month. The corresponding degree 
of the polynomial fit is chosen sufficiently high so as to maintain accuracy com- 
parable with that of the JPL data. The improved efficiency is provided in terms 
of reducing the computational time required for the various calculations and 
transformations which may be required for use by GTDS as described in Section 
2 . 1 . 
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SECTION 2 


GTDS SLP EPHEMERIS FILE GENERATION 


2.1 PURPOSE 

A SLP ephemeris file in GTDS contains data in the form of Chebyshev polynomial 
coefficients. Use of these data (as opposed to the direct use of JPL data) re- 
duces the computational time required for the following calculations and 
transformations : 

a. The calculation of the positions and velocities of planetary bodies. 

b. The calculation of the equation of the equinox and the A.l to UT1 time 
conversion for the Greenwich Hour Angle. 

c. The transformation from (or to) the mean equator and equinox of 
1950.0 to (or from) a true equator and equinox system. 

d. The transformation from selenocentric to selenographic coordinates. 

The above calculations are required for the non-analytic orbit generators and 
for the calculation of observations, while the transformations may be achieved 
for the same purposes if the need exists. The calculations and transformations 
are achieved by evaluating the ephemeris polynomials. The expression for the 
polynomials is given by: 


X(t) = £ 
i = 1 


t- t. 


i - l 


m 


186,400, 


( 1 ) 


where 

X (t) = the desired vector, equation, or matrix to be solved (matrices are 
used in transformation computations) 

kj = the ephemeris polynomial coefficients 

t = the time in seconds from the beginning of the ephemeris year 
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t m = the time in seconds from the beginning of the ephemeris year to 
the midpoint of the curve-fit interval 

n = degree of the polynomial plus one 

Details pertinent to the calculations and transformations given by b., c., and d. may 
be found in the document GODDARD TRAJECTORY DETERMINATION SUB- 
SYSTEM MATHEMATICAL SPECIFICATIONS (March, 1972). 


2.2 DISCUSSION OF GENERATION PROCEDURE 

The generation of a SLP ephemeris file in GTDS is an operation performed by 
the Data Management Program which primarily is responsible for generating 
working files of data to be used immediately by other programs in the system. 
The Data Management Program may also operate as a stand-alone program 
whereby files may be created for future use. The generation of a SLP ephemeris 
file is initiated via the Data Management Program’s testing of switches which 
determines the data management operations to be performed. A determination 
is then made as to which of the following three operations relating to SLP ephem- 
eris file generation is to be performed: 

1. Generate a SLP ephemeris tape file from a JPL ephemeris tape. 

2. Generate a SLP ephemeris disk file from a JPL ephemeris tape. 

3. Generate a SLP ephemeris disk file from a SLP ephemeris tape. 

The creation of any one of these files is initiated through the use of the following 
SLP option parameters which are stored in labeled common as a result of being 
specified via card input: 

1. The starting date of the file. 

2. The number of data records to be created. 

3. The bodies represented in the file. 

4. The degree of the coordinate transformation curve-fits. 

5. The degree of the curve-fit of the ephemeris of the fast moving body. 

6. The degree of the curve-fit of the ephemeris of the slow moving bodies. 
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7. 


The number of days represented by the curve-fits. 

8. The coordinate system reference, i.e., 1950.0 or true-of-date. 

The fast moving body is the non-central body whose rate of change in position 
around the central body is greater than the other non-central bodies, i.e. , slow 
moving bodies, requiring a greater degree for the polynomial coefficients for 
the position coordinates. 

In generation operations 1 and 2, a direct application of the curve- fitting tech- 
nique in generating a SLP ephemeris file is made. The operations begin with 
the writing of the SLP ephemeris tape or disk file header record reflecting the 
values of the SLP option parameters that are stored in labeled common. 


Curve- Fitting with Chebyshev Polynomials 

An initial call to the SLP ephemeris generation control subroutine SLPEPH is 
made to initialize quantities necessary in the curve-fitting process and to com- 
pute the time arrays for a later call to SLPEPH. Subsequently, calls are made 
to subroutine CHEBY which houses or treats the following basic formulation of 
the Chebyshev polynomial curve-fitting technique: 

Chebyshev polynomials of degree m are defined by 

T m (x) = cos (m arc cos x), (2) 


for 


x in [-1, 1 ] . 


The first and second Chebyshev polynomials are 

T 0 (x) = 1 and Tj (x) = x. (3) 

Other Chebyshev polynomials are readily obtained by using the recurrence 
relation 

T m+ 1 ( x > = 2 x T m (x) - T m -i (x) (4) 

Using equations (3) and (4) a Chebyshev series of degree m could be constructed 
in the form: 


Y m (x) = c 0 T 0 (x) + Cj T, (x) + . . . + c m T m (x) 


(5) 


2-3 


Suppose that n points (t lf yj), (t n , y n ) lying in some interval [a, b] are 
given and we wish to find a curve approximating these points which is of the form 
of equation (5). The points must first be changed to the set of points Xj , (i = 1, 

2, . . . , n), for which the Chebyshev polynomials are defined by using the 
transformation 


2tj - (a + b) 


b - a 


( 6 ) 


Applying the least squares criterion at this point maintains that the expression 

2 


n 

s - £ 

i = 1 


m 


Yi" £ c j T j < x i> 
j = 0 


(7) 


should be a minimum. To find the minimum, partial derivatives are taken with 
respect to the cj , obtaining the system of simultaneous equations 

4- = 0, j = 0, 1 Ill (1 

a 'j 

which can be written in the following matrix notation: 

[T] [C] = [P], 


where 


2 Tj <x,) 

2 T 0 ( Xi ) Tj (xj) . 

• ^ T 0 (Xj) T m (Xj) 

2 T, (Xj) T„ (x,) 

2 Tj (Xj) 

• ^ T i (V T m 


[T] = 


s T m (x,) T 0 (X|) 


S T m ( Xj ) T, (x,) r . 2 Tj, (x,) 


n 

where 2 signifies , 
i = 1 
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~ c o 


i 

M 

H 

o 

V 

JX 

1 


c i 


2 Yi Tj (xj) 

[C] = 

• 

and [P] = 

- 


_ Cm _ 


_ 2 Y, T m (xj) 


The need to solve simultaneous equations can be eliminated by taking advantage 
of the following property of the orthogonal Chebyshev polynomials: 


If k and fi are non-negative integers and 

are not both zero then, 


( 

0, fork ^ £ 

(10) 

m + 1 1 



S T k (X,) T, (J,) = l 

(m + l)/2, for k = £ -=h 0 

(ID 

m + 1, for k = £ = 0 

(12) 


where 


(2i — 1) 7T 


(13) 


The use of these nodal points yield the classical Chebyshev polynomials, and we 
have as opposed to relation (7), 

m + 1 

S- E 

i = 1 


Vi- 


nt 

E 

j = o 


C; T; 


j (Xj) 


(14) 


Taking partial derivatives with respect to the Cj in this instance results in the 
diagonalization of [T ] . Thus we have 


[D] [C] = [R] or 


(15) 



’ 2 Tq (X 4 ) 0 

~ 

0 


V 


2 yj T 0 ( Xj ) 

0 ST^Xj).. 

. 0 


C 1 

— 

SyiTj (xj) 

0 0 

2T m (*i)_ 


_ c m 


2 9i T m <*i> 


m + 1 

where here, 2 signifies ^ . 

i = 1 

Thus, we can now readily compute a value for each coefficient Cj , 

m + 1 

£ Pi T j (*i) 

C j =± ^ T~l ’ j 30 ’ 1 ’--’" 1 ( 16 ) 

£ T/txj) 


Using relations (11) and (12), we obtain 


m + 1 


C; = 


J m + 1 


£ V\ Tj (Xj), j -h 0, 


(17) 


i = 1 


and 


C; = 


1 


m + 1 


m+ , E yiW’ j = °- 

1 = 1 


(18) 
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Having determined the coefficients cj , the Chebyshev expansion of degree m for 
Y m (5c) has been completely determined: 

m 

Y m (x) = £ c j T j^>' (19) 

j = 0 

At this point, Y m (5c) has been expressed as a Chebyshev series in the interval 
(-1, 1). The Chebyshev series may be converted to its equivalent power series 
in (-1, 1), and this series in turn may be scaled to the interval (a, b) by using 
the transformation 


(b - a) x b + a 

t = 2 - — + 

2 2 

The resulting power series may be written in the form 

m + 1 

Y m (t)= £ v 4 - 1 , 

i = 1 


( 20 ) 


( 21 ) 


where the A,'s are directly obtained from the Cj’s as follows: 

A! = 2 a ik , i = 1, 2, . . . , m + 1 
i 

k = j + 1 = 1,2,... (22) 

i + 2 (k- 1) < m+1 

where 

a ik = (- 1 ) k + 1 c 2k-i> k = j+l,2,... (23) 

2k - 1 < m + 1 


a^ = 2 1 2 Cj, i = 2, 3, . . . , m + 1 


(24) 



and 


= c i+ 2 (k— 1 ) ^ 2a i-l , k ~ a i, k-1 
i = 2, 3, . . . 

k = j + 1 = 2, 3, . . . 

i + 2 (k- 1) < m + 1. 


(25) 


A pplication in Creating a SLP File 
Given: 

m — degree of polynomial to be fit 

y — array of JPL ephemeris data points (m + 1 components of nutation, 
position, or velocity) 


x — array of times in the interval (-1, 1) that are used to get the data 
points 

[ These times are a function of m and are computed using equation (13). 

cos(2i - 1) n 


X ‘ 2 (m + 1) 


i = 1, 2, . . . , m + 1 


] 

(13) 


The coefficients (c 0 , c lt ..., c m ) are computed first using equations (17) and 

( 18 ). 


C j m + 1 


m + 1 

E y-x T o <*i>’ J # 

i = 1 


0 ; 


(17) 


and 


_ 2 
C i m + 1 


m + 1 


E y; T j (x i)’ j = °- 

i = 1 


(18) 
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The expanded Chebyshev series 


E c j T jW 

! n 


is then converted to its equivalent power series, which in turn is scaled from 
the interval (-1, 1) to the interval (a, b) taking the form 

m + 1 

Y m (0 = E \ «i-l - <2<S) 

i = 1 


where 

t c is the time at the center of the fit 
tj_j is within the limits of the fit 

| b - a | = length of fit and a and b are the minimum and maximum values 
that (t j_j - t c ) can assume 

Aj's are the various polynomial coefficients resulting from the calls to 
CHEBY 

Since I t - t c I < 1/2 length of the fit, 

- 1 12 length of fit < t - t c < 1/2 length of fit 
a b 

These limits are constant during the entire time span. 

The remaining subroutines and functions that are used in the generation proce- 
dure are given in Appendix B, and the overall flow of the procedure is provided 
by Appendix C. 
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SECTION 3 


PROCEDURE FOR USING THE SLP FILE 
3.1 INTRODUCTION 

The SLP file contains the polynomial coefficients for the position and velocity 
coordinates of the fast moving body and the position coordinates for one to seven 
slow moving bodies. The header record of the file contains: 

1. Day of the first record on the ephemeris file. 

2. The year of the ephemeris file. 

3. The number of records (days) on the ephemeris file. 

4. The bodies represented in the file; 

a. central body 

b. fast moving body 

c. one to seven slow moving bodies. 

5. The degrees of the polynomials for 

a. rotation matrix (matrices) 

b. fast moving body position 

c. fast moving body velocity 

d. slow moving body position. 

6. The number of days per curve-fit. 

7. The coordinate reference indicator for the SLP ephemeris (1950.0 or 
true-of-date coordinate system). 

8. The number of bodies to be processed for the SLP ephemeris. 

The remaining records on the SLP file contain: 

1. The time in seconds from the start of the year of the SLP file to the 
midpoint of the day of the year of the SLP file. 
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2. The polynomial coefficients for the position coordinates of the fast 
moving body. 

3. The polynomial coefficients for the velocity coordinates of the fast 
moving body. 

4. The polynomial coefficients for the position coordinates of the slow 
moving body. 

5. The polynomial coefficients for the selenocentric to selenographic 
transformation matrix. 

6. The polynomial coefficients for the 1950.0 to true-of-date transforma- 
tion matrix. 

7. The polynomial coefficients for DELTA H (calculation b. of Section 

2 . 1 ). 

8. The day of the record. 

The SLP ephemeris file is created by the program CREATE3. 

The program will read the JPL DE-19 direct access ephemeris tape, compute 
the polynomial coefficients and write them onto the SLP ephemeris file. 


3.2 JCL REQUIREMENTS 

1. // EXEC PGM=CREATE3, REGION=250K 

This step will retrieve the program to generate the SLP file. 

2. //STEPLIB DD DSN = TESTLOD, UNIT=2321, 

// VOL = SER = RTAC, DISP = SHR 

This step locates the program within the system. 

3. //FT05F001 DD * 

File 5 will contain the input data cards describing the generation of 
the SLP file. 

4. //FT06F001 DD SYSOUT = A 

File 6 is used by the system for printer output. 
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5. //FT14F001 DD DSN = XXX, UNIT = DISK, 

// DISP=(NEW, PASS), SPACE = (CYL, (4, 1» 

File 14 is used only when a SLP file is to be created on disk. 

6. //FT33F001 DD UNIT = 9TRACK, 

// DCB = (RECFM=VS, BLKSIZE = 3460), LABEL=(1,BLP) 

File 33 is used only when a SLP file is to be created on tape; 

7. //FT34F001 DD UNIT = 2400-9, 

// LABEL = (, BLPj , IN), VOL=SER = ABC, 

// DISP = OLD, DCB = (RECFM=VS, LRECL = 8304, 

// BLKSIZE = 8308, DEN = 2) 

File 34 describes the input DE-19 JPL ephemeris tape. 

8. //FT60F001 DD DSN = GTDS .ODS. ACCOUNT, 

// DISP = SHR, UNIT = 2314, VOL=SEft=GlUSRl 

File 60 contains the accounting information for the GTDS system. 


3.3 INPUT REQUIREMENTS 

To generate the SLP file a DE-19 JPL ephemeris tape and input data cards are 
necessary. The JPL tape is described by file 34 and the input data cards by 
file 5. The data is input through NAMELIST. 

The format of the input cards is: 

SLPEPHEM 

&SLPDAT 

Option cards as described on the next page 
&END 
FIN 
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DATA FOR SLP EPHEMERIS FILES 

CARD TYPE FORMAT 

I. DATA &SLPDAT NAMELIST 

IOPER Operation to be performed 

= 3 create SLP disk file from SLP tape 
= 4 create SLP tape from JPL tape 
= 6 create SLP disk from JPL tape 

SLPYMD UTC Calendar date of start of file 

(YYMMDD. ) 

ISPAN Number of data records to be created 

{Default = 1, maximum of 11 on disk file) 

NBEPM(I),I=1, 9 Body numbers of bodies represented in the 

file (default = Earth, Moon, Sun, NDEGRE(1)=1 
NBEPM (2) = 2, NBEPM (3)=3) 

1=1 central body 

=2 fast body 

=3,9 slow bodies 

The valid NBEPM settings are: 

=1 for Earth =6 for Saturn =11 for Venus 
=2 for Moon =7 for Uranus 
=3 for Sun =8 for Neptune 
=4 for Mars =9 for Pluto 
=5 for Jupiter =10 for Mercury 

NDEGRE(I),I=1,4 degrees of curve-fits (Defaults; 

NDEGRE (1)=4, NDEGRE(2)=8, NDEGRE(3)=4), NDEGRE(4)=8 

1=1 degree of coordinate transformation matrices 

fit (also used for GHA correction term for 
A.l to UT1) 

1=2 degree of position fit for NBEPM(2) (fast body) 
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1=3 


degree of velocity fit for NBEPM(2) (fast body) 


1=4 degree of position fit for NBEPM(3), ...» 

NBEPM(9) (slow bodies) 

NCFDAY Number of days represented by each curve- 

fit (Default=l) 

ISLP50 Reference of SLP coordinate system 

(Default=l) 

=1 for 1950.0 
=2 for true-of-date 

&END 

Note: The JCL requirements for the SLP data sets are: 

IOPER FORTRAN data sets used 

3 FT14F001, FT33F001 

4 FT33F001, FT34F001 

6 FT14F001, FT34F001 

where: 

FT14F001 Defines the SLP ephemeris disk file 

FT33F001 Defines the SLP ephemeris tape file 

FT34F001 Defines the JPL ephemeris tape file 

An example of the use of CREATE3 in generating a SLP ephemeris disk file is 

given on the next page. 
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//CMIBTSL3 JOB | QN«? I 4 I 3 1 IM • T *G 0 ->4 35 , 'VJI 9*3 I , CCC 

✓ ✓ exec PGM S rp CA T 63 «PRCi !ON= 2tCVK 

//STEPLIU DO OSMsTfcSTLUC.UMIT = 232l#VOL=SfcR = RTAC*OI SP=SHR 
✓/FTG5FC9I 00 * 

✓ /FTf-6Frot OO SY SOU T = A 

//FT13FPC 1 DO Oi>N = GT(>S.CDS.enni)PMSGt UNI T = 2JI4 *D I SP =SHR * V 0L=S5 R=G 1 USR 1 
//FT1AFCC 1 90 USN=CNF*SLP«UNI T=DISK*DI SP = ( Nt W# P ASS > *SPACE = (TRK t <lC*l>> 
/✓FT 3 AFP 01 00 UN I T= £40 0-9 .LAf;FL=< »9LP» .IN) . VOL= SER=33 9 1 2D *D ISPsOLO* 

// DC0=( PECFM=VS.LP E CL = R 33 A ,ULK3 I L\L =Q3C9* 35 N=2) 

//FToOFOCl DU DSN=G TCS*OOS*ACCJUNT *D I SP=SMR • UNI T =23 1 4 • VOL = SER=GI USR 1 
//NUCLEUS DO O IGP^SHO.UNI T=S YSDA . VOL=PEF=S YS I . S VCLI 0 

IEF236I ALLOC • F CR CM I (3TSL 3 I 

IEF237I 2EJ/0 ALLOCATED TO STEPLI'J 

I EF 2371 231 ALLOCATED TO FT05TODI 

IEF237J 33A ALLOCATED TO FTCEFOOl 

1 EF 23 71 IA2 ALLOCATED TO FT13F001 

IEF237I J3A ALLOCATED TO FTI4F091 .. .. . - 

IEP237I CC 2 ALLOCATED TO FT3AF>jJ 

IEF237I 14 2 ALLOCATED TO FT69F00I 

IEF237I ICC ALLOCATED TO NUCLEUS 
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— GTOS PUN NUMBER IS ie?S 
STATING ADDRESS OF MAIN IS 2376902 

— COMPUTER IDENTIFICATION l S Gl 

- - JOB IDENT IFICATICN IS CMIUTSL3 

CURRENT TIME IS P3C0 3* 19 72 A I- A* 2A GUA Y , 47 • 9 SC 00 A3 0 
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-CSLPDAT ' 

I SPANS . IltNBF.PMs • la 2. * 3. 0* r • °* 0 

- NOE ORE* --- 8* 6* 3. 8*NCPDAYa 1 • ISLP5C* ' 1 . IOPER« ® * SLP VMD» 

670523* COOCO^COOO •NDSLf s 3*I0AYl= 1A3*DJ= 24 306 33*500 437654 • I YE AR = 1967 

-fcEND “ 




SECTION 4 


THE USE OF THE SLP FILE 
BY THE STAND-ALONE SUBROUTINE SUNRD 


Once the SLP ephemeris file is created, the stand-alone subroutine SUNRD will 
read the polynomial coefficients for the position coordinates of the bodies in the 
SLP file and compute the position vectors for the bodies requested through the 
calling sequence of SUNRD. 


In put Requirements 

SUNRD requires an SLP file containing the bodies for which position vectors are 
requested, the Julian date, the number of the bodies requested, and the real 
function DJUL to execute. 


SLP Ephemeris File 

The SLP file to be read by SUNRD is described by file 14 as follows: 


//GO.FT14F001 DD DSN = XXXX, UNIT = DISK, 
// DISP=SHR , V OL=SER=XXXX 


If the SLP file was temporarily created for the particular run, the DISP param- 
eter is changed to DISP = (NEW, DELETE). 

When SUNRD reads the SLP ephemeris file, it stores the parameters from the 
file into labelled COMMON, from SUNRD stores the header recorder into 
COMMON /SLPOPT/ and the remaining records in COMMON/SLPREC/. The 
only variables used by SUNRD from the SLP file are: 


a. IDAY1 

b. IYEAR 

c. NBEPM 

d. NDEGRE(2) 


The day of the first record on the ephemeris file 

The year of the start of the ephemeris file 

The central body, fast and slow moving bodies for the 
polynomial coefficients 

The degree of the fast moving body’s position 
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e. NDEGRE(4) 

f. NCFDAY 

g. NBSLP 

h. TSEC 

i. PPOLY1 

j. PPOLY2 

k. IDAY 


The degree of the slow moving bodies’ position 

The number of days per curve-fit 

The number of bodies processed for the SLP file. 

The time in seconds from the start of the year re- 
quested to the midpoint of the day 

The polynomial coefficients for the position coordinates 
of the fast moving body 

The polynomial coefficients for the position coordinates 
of the slow moving bodies 

The day of this particular record. 


THE CALLING SEQUENCE 

CALL SUNRD (DATE , NB, POSVEC , IERR) 

INPUT 

FORTRAN 

NAME DIMENSION DESCRIPTION 


DATE 1 Julian date for which the position vectors are 

to be computed 

NB 9 The number of the bodies requested. 

NB(1) = central body 
NB(2) = fast moving body 
NB(3)-NB(9) = slow moving bodies 
The bodies are represented by the following 


numbers: 


1 = earth 

7 = uranus 

2 = moon 

8 = neptune 

3 = sun 

9 = pluto 

4 = mars 

10 = mercury 

5 = jupiter 

6 = saturn 

11 = venus 
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OUTPUT 


FORTRAN 

NAME DIME NSION DESCRIPTION 

Position vectors for the bodies requested 

Error indicator 

0 = no error 

1 = I/O error 

2 = error in specifying the central or fast 
body 

3 = non-central body not in the SLP file 

When an error is encountered a message will be written and the corresponding 
error number will be transmitted through the IERR parameter. If no error is 
encountered, the number 0 is transmitted. 


POSVEC 3 » 8 

IERR 1 


DJUL 

The real function DJUL is used to compute the modified* Julian date of a given 
Gregorian date after 1950.0. 

The following seven pages consist of a listing and flowchart of SUNRD and a 
listing of DJUL. 


♦Modified Julian Date = Julian Date - 2430000. 
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LEVEL £0 • 1 I A UG 71) (.5/160 FORTRAN H - 

.... COMP ILf.ff OPTIONS - NAME - M A I N • C PT = 0 i *L I NECNT =5« , 5 I Z E=OOOOK . - — — 

SOURCE. F.eCOIC.NOLIST.NCOECF .LCAC.MAF.NOEOIT. ID. XRfc'F' 

1SN OOC 2 - SUBROUTINE SUNMO ( DA TE . M). PCSVE C . 1 1 P R t 

C*4**44*4******4***4W.Y*444***4*> !; **4W**W(W**'*»F4F*A4F*F* , ' , F*****F*M****** 

c 

- - - C Vc PSION OF Ort/Ol/72 — — - 

C VERSITN OF 11/9/71 

C - VERSION CF 1 H! /I / 72 ■ - — 

c POP TRAN M SUBROUTINE FCP I EM S/360 

c PUPPC SF - - 

C TO COMPUTE THE POSITION VECTORS 

C CALLING SEOUENCE - — - - 

C CALL SUNROICATE .NO.PCSVEC. IERR) 

C CALLING SEQUENCE 

. C — ' - DATE - - JULIAN UNIVERSAL CATE 

C NO THE BODIES RE GUEST E C 

C - NE (1 ) - CENTRAL 60CV - 

C NG (2 ) - THE FAST MOVING etJOV 

C NHI3-FKI — THE SLOW MOVING-GOO IES 

C OUTPUT 

C CALLING SEQUENCE - 

c posvec - the position vectors fcr tfo rodies specified 

c ■ -- ... - - - - IN THE NE ARRAY — 

- c 

o - -IFRR --ERRCR INDICATOR 

C 0 - NC 0RRCP 

C _ 2 - ERROR IN SPECIFY ING-TNE CENTRAL-OF FAST -- 

c bcoy 

..... .... c 3 - NON— CENTRAL ECDV NOT IN TNG SLP-FIL6 t 

— C COMMON NLCCK INFORMATION USED INTERNALLY 

C COMMON /SLPCPT/ 

..... ... c IDAY1 - CAY (IF FIRST RECCRO CN EPFEMERIS F ILE 

C - I YEAR - YEAR CF TFE EPFEMERIS FILE 

— ... C - NBEPM - BODIES FCR PCLYNCMIAL COEFFICIENTS — 

C NQEFMdl - CENTRAL ECCY 

{ - NdEF M( 2 ) - FAST MOVING POCY 

C NJEFMC3FK) - SLEW MOVING RODIES 

- C NOEGRE - DEGREE CF POLYNOMIALS — 

C ' NOEGRE 12) - FAST MOVING HOOY POSITION DECREE 

- — C NOEGREI&) - SLOW MOVING eOC I E S -POS I T:iON DEGREE- 

C NCFOAY - NUMBER OF CAYS PER CURVE FIT 

c NBSLP - NUMBER OF ECCIES CN = ILE 

1 C COMMCN /SLPRCC/ 

- c TSEC TIME IN SECCNCS FRCM START OF-THIS YEAR 

C TO MIDPOINT CF T F IS RF.CORC TIME INTERVAL 

— -c PP0LY1 - THE PCLYNCMIAL COEFFICIENTS FOR THE FAST ; — 

C MC V I NG BCOY 

C - - - PPCLY2 - FCLYNOMIAL C CEFF 1 C IENTS -FOR -TFE — SLOWF 

C • MOVING BODIES 

.... c ICAY - BEGINNING CAY CF THIS RECORC - 

C 

- — C*#*»*»**.«<***4******«****4»**«4#*«4i***-*4***ii*Y«»****4*******-*«*»*****4» • 

c 


DATE 72.339/17.43.03- 
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rfA. 

I 

Ol 


I SN 0003 
ISN 0CC4 

- ISN 0008 

ISN C0C6 


ISN C0C7 

isn coce 

ISN 0009 
ISN 0010 
- ISN coil - 

c 

c 

ISN CO 1 2 

c 

- c • 

c 

-ISN 0014 


ISN 001* 
ISN C01« 


ISN C01 7 
ISN 0018 


ISN COl S 

isn co 2 o 

ISN 002 1 
ISN 0023 


ISN CO 2 4 


c 

c 

c 

c 

c 

c 



c 

c 

c 


ISN C025 
ISN 00 20 — 

ISN CC27 
ISN 0028 
ISN C 02 9 
ISN 0030 
ISN CO 2 5 

ISN 0033 

ISN CO 25 
ISN C03C 
ISN CO 3 7 


IMPLICIT REALMS (A-H.O-Z) .• • : - r 

LOGICAL FIRST 

CUMMCN/SLPCPT/ OJ vICAVt • IYCA R ■ tSPAN — - NB EPM ( 9 I - ■». — — • 

* NCEGREIAI • NCF DAY • ISLP 50 .NHSLP 

COMMCN/SLPREC/ TSEC . FPCL Y 1 < 3 . 20 ) - . VPOLY 1 < -3.-12 > — - 

1 FFOLY2I3.13.7I , AFOLYI 3, 3, 10 ) • 

2 CF0LY(3.3.10I .PCELM 10), ICAV - — 

DIMENSION EFT I N.E <20 ) • PB <3,81 ,FCSVEC< 2. 8 ) .NB( 9 ) t 

1 - — ■ — 1SLB<9) - ' — ; 

DATA I OF RE V/ -99999 99 9/ »F IRST/,T RUE ,/ .NWSLP/14/ 

DATA TOL / 1.00- 30/ - ••••• — - 

DEFINE FILE 14< 100 ,1 1 27 , U, I AV ) 

l ERR = 0 -- : . 

IS THIS THE FIRST ENTRY-INTC-SURIRD KC — . — — 

•• 

IF (.NOT. FIRST) GC TO 2 — r 

RE AO EPHEMERIS HEADER READ AND STCRE IN-COMMON/SLPOPTA-. 

READ IN.SLFU .ERR = eOO) 1C AVI, IYEAR, ISPAN.-- NBEPP . 

•NOFOPE, NCFOAV. ISLP50 . fcHSLP 


FIRST = .FALSE. 

Y = OFLCAT (IYEAR - (IYEAR/100) « 100 ) 


CONSTANT At -UTC OFFSET FCR 1987 


A1UTC = 5 • 662 5 2 ECO 


D J.LLF = 0 J LL <Y.l.CO.l.DO.O.DO,O.CO,O.CO)N’4 13 000. ODJ 


IS CAY IN SANE INTERVAL AS PREVIOUS CAY 


2 lx(OATC -OJLLE)«UC400.000 ♦ AltTC 
IDAYR = IOINTIT / 46 A- 00.000) ♦ I 

IF (IOAYR * Gf: . IOPFEV .ANO. ICAYR ,LT . ( I RPREV 4NCF C A Y ) ) GO TO 3 N, 

IREC = < ( I O A YR - ICAYl) / NCFCAY-) ♦ 2 

READ EPHE PERIS COEFFICIENTS — — - 

READ (N, SLP # I RFC .ERR = -30 3) TSEC . F F CLY I . VPOLY l , PP0LY2 , APOL Y rCPOL Y . --- 

1 PDcLH . ID AY _ 


IOPRE V = I O A Y 

3 EPTIMEII) = < ( T - TSEC) 86400. OCO) - — - 

VAN = maxO (N0EGRG(2) .N06GRE<4)) 

DM 4 I = 2 .FAX - - - 

CP T I ME II) = EPTINE(l) * EPTINF(I-l) 

IF<OARS<EPTIME< l )).L T. TOL) -ERTIMSI I ) = O.OCO 

4 CONTINUE 

IF <NB< I ) .EO.O.OB.NB<2) .EG.O) -GC TO 930- 

N3CEM =99999999 

NAP =0 • - —A— - ~ — 

I = c 
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ISN CO 3 6 K = 0 

I SN C 0 3 9 t C C I = I ♦ 1 

ISN coil) IHd.ro. ID) GO TC 20 ) — ■ ■ 

ISN CC«2 I'INBU I.EO.O) GC TC 20 ) 

ISN C044 NAR = NAB ♦ I 

ISN C04S 00 ISO J=I .NOSLP 

ISN CC46 IF(N6(I> .EC.NBEFM J) » GO TO 17S 

ISN C048 ISC CONTINUE 

ISN C C 4 9 GO TC '30 • - 

ISN COSO 175 IF(J.FQ.l) GO TC 1 HO 

ISN COS 2 - K » K I - - 

ISN COS 3 ISliB(K) = J - 1 

ISN C OS 4 - - GO TC 100 — - - - 

ISN COSS ISC N0CENT = 1-1 

ISN- CCS*- — GO TC 100 - 

ISN CC57 200 CONTINUE 

ISN cose-- NFOOG = NDEGRE ( 2 ) ♦ 1 - 

ISN CC5S NSBCG = NOEGRE < 4 ) ♦ 1 

C EVALUATE POSITION POLYNOMIALS 

ISN 0060 21= DO 22C Jj=t,K 

isn ooe i j = isls(jj) 

ISN CCS 2 IF(J.EO.l) GO TC 200 

ISN COS 4 JH1 = J -1 

ISN CC65 Bfld.JJ) = PPCLY2U .1 ,JM1) 

ISN 0066- 66(2. JJ) = PPCL72 (2 .1 , JW1 ) 

ISN CC 6 7 06(3. JJ) = PPCL Y20 .1 . J»1 > 

ISN CC6S DO 220 I=2.NSE)0C 

ISN CC69 FIQ(I.IJ) = mid .JJ) PPCLY2 <1 . I . J»t ) - =PT IMF! I- 

ISN C07 0 60(2. JJ) = H0C2.JJI ♦ PPCLY2 (7 . I . JN1 ) * SPT IMS) I- 

ISN CC7 l 6HC3.JJ) = jl(3»JJ> N PPCLY2 ( 3 , I . JMl I * 5PTIMEII- 

ISN 0072— GO TO 220 - • 

ISN C C 7 3 23 0 H.l(l.JJ) = PPCLYld.t) 

ISN CC74 63(2. JJ) = PPCL1K2.1) - - 

ISN C07S 66(3. JJ) = PPCL7K3.il 

ISN C C 7 € DO 2«0 I=2.N=Ui)G 

ISN CC 77 06(1 .JJ) = 36(1 .JJ) ♦ PPOLYl(l.I) 4 EPTIHEII-t) 

ISN CO 76 - - 66(2. JJ) = 66(2 .JJ) ♦ PPCLY1 (2.1) 4 5FT IMc( I— 1 ) - 

ISN C C 79 60(3. JJ) = 63(3. JJ) ♦ PFCLYl (3.1) 4 EPTIMe(I-l) 

ISN cceo 24C CC N T IN LE 

ISNCCet 22 'C CONTINUE 

ISN CC6 2 K = 1 - 

ISN CC63 NAP = NAR - 1 

ISN ■ CC6 4 — - IF (NBCENT.EQ.O) K = ■> - - - - 

C 

C OUTPUT POSITION VECTORS 

C 

ISN CCS 6 DO 260 J = I .MAR 

ISN C 0 E 7 IF{ J-N6CENTI270 .280.290 

isn ccee 27 c pusvec d . j> = ned . J4-1 » - eoci ,1 1 — - - 

ISN CC“9 PUSVFCI2.JI = je(2,J.I) - RB(?.1> 

ISN C090 ~ PO S VEC ( 2 * J I = nen.JM) - 0B(3,I) - -- 

ISN CCS 1 GO TO 2 6 0 

ISN CCS 2 - 2 = C PO S VEC ( 1 . J) -= -5E (1 .1 h - - 

ISN CCS 3 POSVECI2.J) = -EB(2.1) 
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O0CO02S7 
CCCCC30 7- 
00 00 0 3C 3 
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ISN 0054- 
ISN CCS 5 


ISN- CO 96 

ISN CCS? 

-590 

ISN - COS 4 

- •- - -■ 

ISN CC99 

sec 

ISN 0100 

. 

ISN 0101--— 

c 

- - £00 

ISN C 1 02 


I SN 01 03 

1400 

ISN 0 I C 4 


[sn cioe — 

520 

ISN ’cics 


ISN CIO?- - 

ICOO 

ISN 0 0 8 


ISN 009 - 

53C 

ISN OlO 
ISN bill— 

1CC1 

ISN 0 12 

595 

■ *;• 

— c - • 

ISN 0113 



-K*t 8!1 ,I>~~ 
K*u 0 ! 2 . 1 ) 
K*H0(3tll- 


posvf.c! 3 -ee<3 t t > 

GO TO 2E0 

posvec! i-.j>-=-bqii . J» ■ 

p(>svfC(2tJ) = ae<2.j) • 

POSVEC<3,J) =-.je<3tJ> • 

CONTINUE 

GO TO ««« • - ; 

HP 1 TE ( 6 • 140M — 

ItRR =1 

FORMAT! *CI/C ERFCF'l \ — 

GO TO 59 5 

WKI TE16. 10901 — — : - - 

IEHR =2 

FORMAT! *0 -ERROR-IN SPEC IFV I NG- TNE -CENTRAe- -AND/OR FAST- 0OCV* 1 
GO TO 599 

»RI TE (£. 190 I - - - • ; — - 

IERR =3 

FORMA Tl • C NCN-CENTFAL EdOIES NCT IN FIUE.*) - 

RETURN 

*«************* 4 »******»- - — 

END 


CCC00315 


- -GQ00031G 
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(7) RETURN J 
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LEVEL 


ISN 


- ISN 
ISN 

- ISN 
ISN 
ISN 

ISN 

-ISN 
ISN 
- ISN 


£0.1 (ALG 711 


CS/360 FORTRAN H 


-DATE -■-72»3-39/l***»-.30- 


COO 2 


COMPILER OPTIONS - NAVE* MAI A * OPT*0 i »LINECNT-58#51ZE*OOOOK. — 

SOUPCE .EBCDIC. NOLIST. NCOSCF ,LCAC. MAP. NOEC tT . ID, XREF 

REAL FUNCTION 0JUL*8 < Y.XM .O.FR .TM, SSCli DJIA- t 

- . . - . ... -OJUL 2 -- 

....... ...i.... ........ ........ . ..OJLL 3 

VERSION CF 2/30/71 OJtX. 5 

FORTRAN FLNCTICN SUBPROGRAM OJUL 6 


C 

c 

c - 
c 

--C- 

. c 

- c 
c 

- c 
c 

- c - 
c 
c 
c 
c 
c 

- c 

c 

■c 

c 

c 

c 

... c ■ 

c 

c 

c 

c 

-~e*- 
• c 
c 
c 
c 
c 
c 


PLRPOSE 


TO CCMPUTC THE NOOIFIEC JULIAN CATE OF A GIVEN GREGORIAN 
CATE AFTER 19SO. 9 (JULIAN LATE - 2430000-* 


REAL FUNCTION 0 JUL< V .X A • 0 • HR . TM . S E C > 


OJUL 7 

OJUL- 0 

OJUL 9 

OJUL 10 

OJLL 11 
-OJL-L 12- 


INPLT 

X - YEAR — • ••• - • 

X M - MGNTH 

•O'- CAY ; 

HR - HOLRS 

T» - MINLTES •— 

SEC - SECONDS 

OLTPLT 

- - -OJUL - THE MODIFIED JULIAN CATE 

REMARKS -■ • - 

INFLT OATE MUST OCCURE AFTER 1950.7 

SLPRCLTINES RECL1KEO 

- NCNE •• - 

programmer 

M. MC GA-1FY — - 


1 3. 

- 1 4— 


OJLL 
OJUL — 

OJUL IS 

OJUL 16 

DJUL 17 

OJUL 13 

OJUL 19 

DJUL- — 20 

DJLL 21 

— OJUL 2 2- 

OJUL 23 

- DJLL 2 A 

OJUL 25 

-DJU 26 

OJUL 27 

OJUL- 28 

OJUL 29 


0003 
CC04 
CSC 5 
C0C6 
C 0 C 7 

ooce- 

0009 

0010 
0011 


• ••••••••••••••«•••• •• • •••••••••••••••■ ••••••• •▼♦r* * • • 

IMPLICIT REAL** (A-H.0-2) -■ 

I =YF1900.0D0 

K -0 

1 I (HR - 12.000 * 3603.007 + TM 4 60.0C0 ♦ SEC) / 8SA0C.3C0 
O JUL=K- 3 20 7 SF 1461 * (IffSOO-MJ-JAI/IS 1/4F 357* (-J-2— <-J— 1-4-1/-1 2* If 1/ 1 
**( ( I F450CF ( J- 14 )/12)/130)/A- 24 303 CO 

OJUL = OJLL F FE - - - 

RETURN 

END - - — - - - 


OJUL 30 

DJUL 31 

- OJUL— 32 

DJUL 33 
.-.-..DJUL- — 34- — 
DJUL 35 

DJUL — 36 

DJLL 37 

— OJUL 38- — 

DJUL 39 

— DJLL- - 40 

DJLL 41 

2— 3DJLL- 4-2 

DJLL 43 

OJLL— 44 

DJUL 45 
-OJUL 46- 



The following is an example of a DRIVER for SUNRD. This DRIVER indicates 
the bodies and date for which the position vectors are requested. After each call 
the position vectors are written out by this driver. 
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LEVEL 20. I CALG 51 


03/361 FCRTRAN H 


...... - COMP ILER CPTI ONS - NAME = V AI N « CPT=0 1 »L I NE CNT = 58 • 3 I i E=OOOOK » — 

SOURCE .EECOI C. NOLIST .NCOECN .LC AC . MAP. NOEC IT . 10. : 

ISN 0002 - — -IMPLICIT PFAL «6<A-H.C-JC> 

ISN 0002 UCFINE FILE 1 A < 20 . 1 1 27 ,U . I I » 

ISN CC04 - U 1 MF. NS l C N NSC9I .PCSVECC3.«> 

ISN CCCS DO 4 J = t .4 

ISN CCCS 4 NBCJ)=0 

ISNCCC7 N8 < 1 I = 1 

ISM coce - N3C2I ■> -s — 

ISN CCCS NH ( 3 ) = 0 

ISN 0010 CATE = 2479633.500 

ISN 0011 CALL SUNFDCDATE .Nb.PCSVEC. IEFRI 

ISN 0012 - *MI TF<6.1?.34)FCSVEC 

ISN C01 2 OATE = DATE .SCO 

.—ISN C014 CALL SONRD {DATE .NO .PUSVGC i I E PR I 

ISN 0015 WXI TE<6.1234)FCS\.EC 

ISN COle — OATE = DATE ♦ .SCO - - - - 

ISN 0017 CALL SLNROIOATE.NO.POSVEC.IEr.pl 

ISN 0010 l»RI TFC6.123AIFCSVEC 

ISN 0019 NUIII » I 

• -ISN C020 - ■ — NO < 2 I =■ 2 - - - 

ISN C05 1 Nrt ( 3 I = 4 

ISN C022 NO(a; =11 

ISN 0023 DATE = 2434033.500 

ISN 0024 .... CALL SUNRDCOATE .NG.PGSVEC.IERfil 

ISN 0025 WRITE IF .12341 FQS VcC 

- -ISN C02« -- CATE = GATF ♦ .SCO 

ISN 0027 ‘ C<LL SLNHO (DATE iNB.POSVEC. IEFR) 

ISN 0028 l»r( S T. < 6 , 1 2 34 | PCS VEC 

ISN 0024 DATE = It A TF. «• .ECO 

ISN CO 30 CALL SL-NPD {DATE .NO .PCSVEC. IEBB) 

ISN 0031 WR I TE ( a • ) <’W > FCS VEC 

ISN 003 2 - NS Cl > = A 

ISfC 003 3 NO <21 - 1 

- ISN CO 3 4 NO (3 ) = 11 

ISN 0035 NH ( 4 > = 7 

ISN CC 3C Nr)< 5 ) = 2 

ISN 0037 OATE = 2 4 2 C 6 33 • 50 0 

- ISN 002H •- CALL SUNFDCOATE .NH .FOSVEC. IEFF 1 

ISN C 029 WRI TEC£.1334)FCSVEC 

-|SN COAO DATE = DATE + .EDO 

ISN 004 I ** CALL SUNR 1 ? COA TE • NO .F CSVEC. I ERR I 

ISN C042 WRI TFC6.1234IFCS4EC 

ISN C043 OA TC = OATE 4 .500 

ISN C044-- CALL SUNRD <0 ATE « N(3 .PCSVEC • I E FF I 

ISN‘0045 WKI TEC6,1P34)FCSVEC 

ISN 00*6 123 4 F OPM AT ( *0* . 1 O X , • FCSVEC* m/ ,2X .8 <✓ . 3X . 3 CSX . 01 7 .8 » I |- 

ISN 0047 RETURN 

ISN -CC4e END — 


- c °py. 



SECTION 5 


ACCURACY OF THE CHEBYSHEV REPRESENTATION 


Specifying the earth as the central body, the moon as the fast moving body, and 
the sun as the slow moving body, the accuracy of the Chebyshev representation 
of the lunar ephemeris was examined as a function of the two SLP option param- 
eters specifying: 

1. The degree of the curve-fit of the lunar ephemeris (values of 6 through 
19 were used) 

2. The number of days represented by the curve-fits (values of 1 through 
28 were used) 

Other SLP option parameters were set as shown in parenthesis: 

• The degree of the coordinate transformation curve-fits (4) 

• The degree of the curve-fit of the ephemeris of the slow moving 
body (4) 

• The coordinate system reference (1950.0) 

• The starting date of the file (720101) 

The rms (673 points) errors in position for the polynomial approximations are 
presented in Figure 5-1 and were obtained using 


rms error = 


673 

E (°i- c i) 2 

i = 1 


/673 


Vi 


where 

Oj = the JPL system value (not quite as accurate as 0.2 Km) 

Cj = the GTDS system value 

The rms values represente d by the graphs are the "scaled rms" values. They 
should be multiplied by 1/V 673 « 1/25. 94 to obtain the actual rms values. A 

sample test run is given by Figure 5-2. 
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km 


RMS POSITION ERROR* 
(CHEBYSHEV POLYNOMIALS) 



TOTAL INTERVAL OF CURVE FIT 
(T q = JD 2441317.500486995) 

Figure 5-1. RMS Position Error of Curve Fit as a Function of Total 
Interval and Degree of Polynomial 


Sample lunar position and velocity SLP file accuracies comparable to the JPL 
ephemeris are given in Table 5-1. The results show that the Chebyshev repre- 
sentation is indeed within the accuracy of the JPL ephemeris system for long 
arc lengths of approximately 28 days when the corresponding degree of the 
polynomial fit is chosen sufficiently high. 
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PACE 
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Figure 5-2 



Table 5-1 


Sample Lunar Position*- 1 ) & Velocity* 2 ) SLP File 
Accuracies Comparable to the JPL Ephemeris* 


D 

E 

G 

R 

E 

E 


RMS Error of 
Chebyshev 
Representation 

D 

E 

G 

R 

E 

E 

1 

RMS Error of 
Chebyshev 
Representation 

6 

5 

0.0721 Km* 1 * 

10 

14 

0.1370Km (1) 



0.367 x 10~ 6 Km/sec (2) 



0.140 x 10~ 5 Km/see (2) 

6 

6 

0.2063 

10 

15 

0.2068 



0.120 x 10‘ 5 



0.231 x 10~ 5 

7 

8 

0.1785 

11 

16 

0.1245 



0.119 x 10- 5 



0.121 x lO- 5 

7 

9 

0.3867 

11 

17 

0.2126 



0.259 x lO" 5 



0.218 x 10" 5 

8 

10 

0.1699 

12 

18 

0.1237 



0.134 x 10" 5 



0.117 x 10' 5 

8 

11 

0.4076 

12 

19 

0.2031 



0.322 x 10‘ 5 



0.200 x lO" 5 

9 

12 

0.1740 

13 

21 

0.1898 



0.156 x 10 -5 



0.181 x 10~ 5 

9 

13 

0.2908 

13 

22 

0.3340 



0.279 x 10' 5 



0.327 x 10~ 5 


* Accuracy « 0.2 Km. 

**Arc Length is Expressed in Days. 
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Table 5-1 (Continued) 


D 

E 

G 

R 

E 

E 


RMS Error of 
Chebyshev 
Representation 

• 

D 

E 

G 

R 

E 

E 

1 

RMS Error of 
Chebyshev 
Representation 

14 

22 

0.1194 Km (1) 

17 

27 

0.1556Km (1) 



0.112 x 10" 5 Km/sec (2) 



0.180 x 10 -5 Km/sec (2) 

14 

23 

0.2179 

17 

28 

0.2542 



0.215 x 10 -5 



0.295 x 10“ 5 

15 

24 

0.1509 

18 

28 

0.1191 



0.158 x 10‘ 5 



0.139 x 10' 5 

15 

25 

0.2753 






0.295 x 10' 5 




16 

25 

0.1142 

19 

28 

0.0659 



0.124 x 10' 5 



0.669 x 10' 6 

16 

26 

0.2029 






0.229 x 10“ 5 
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APPENDIX A 


JPL EPHEMERIS SYSTEM IN BRIEF 

The data contained on the JPL DE-19 tapes which collectively span the time in- 
terval from December 30, 1949 to January 1, 2000 may be represented in the 
following tabular form: 


t 

to 

■ 


tf 

r 

r o 

r i 

• • • 

r f 

r 

*0 


• • • 


At p 

© 

< 

At//j 

• • • 

At// f 

Ae 

Ae 0 

Ae i 

• • • 

Ae f 

At// 

D> 

€- 

o 

At//j 

• • • 

AtJ/f 

Ae 

Ae 0 

Ae t 

• • • 

Ae f 

d 2 Q 

d 2 Q 0 

d 2 Q, 

• • o 

d 2 Q f 

d 4 Q 

d 4 Q 0 

d 4 Q, 

• • • 

d 4 Q f 


where 


= fitted integration positions of the planets of the solar system, 
the earth-moon baryeenter, and the earth’s moon. 

= corresponding fitted integration velocities 

= moon nutations in longitude 

= moon nutations in obliquity 

At// and Ae = corresponding nutation rates 

d 2 Q = second modified differences of the above quantities 

d 4 Q = fourth modified differences of the above quantities 

t Q = initial Julian date for which data are provided 

t f = final Julian date for which data are provided 

stepsize for lunar data, At // and Ae = 1/2 day 
stepsize for mercury data = 2 days 

stepsize for all other data = 4 days 


r 

At// 

Ae 
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Positions and velocities are referred to the rectangular equatorial reference 
frame of the mean equator and equinox of 1950.0 = Julian date (JD) 243 3282.423. 

Planetary data are heliocentric and are expressed in astronomical units (AU) 
and AU/day, and lunar data are geocentric and are expressed in units called 
"earth radii" and "earth radii/day". 

The lunar and planetary data were generated by a least-squares fit to source po- 
sitions obtained on the basis of current planetary theories. Velocity coordinates 
for the lunar ephemeris and nutation rates were computed by numerical differ- 
entiation, while planet position and velocity coordinates were obtained as a nu- 
merical integration fit to source positions. The uncertainty in the geocentric 
position of the moon's center of mass is estimated at 150 meters, and the un- 
certainty in the distance is approximately 60 meters. Modified second and fourth 
differences are retained to facilitate the use of Everett's fifth-order interpolation 
formula in calculating intermediate values. The truncation error bounds asso- 
ciated with the interpolation formula are given in Table A-l. Optionally, data 
reduction can be achieved by curve-fitting techniques at the risk of a possible 
loss in accuracy. 


Table A-l 

Bound for Truncation Error 

When Using Fifth-Order Everett Interpolation Formula* 


Body 

Position 

Velocity 

Mercury 

8890.00 AU 

4420. 00 AU/Day 

Venus 

4.73 AU 

0.62 AU/Day 

Earth-Moon Barycenter 

5.19 AU 

2.50 AU/Day 

Mars 

6.74 AU 

5.77 AU/Day 

Jupiter 

6.64 AU 

5.72 AU/Day 

Saturn 

6.64 AU 

5.72 AU/Day 

Uranus 

6.64 AU 

5.72 AU/Day 

Neptune 

6.64 AU 

5.72 AU/Day 

Pluto 

6.64 AU 

5.72 AU/Day 

Moon 

10100.00 Earth Radii 

14500.00 Earth Radii/Day 

At ]> 

0.46 Radii 

1.16 Radii/Day 

Ae 

0.23 Radii 

0.58 Radii/Day 


*A11 entries have been multiplied by 10 12 


The above table is taken from Reference 3. 
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APPENDIX B 


MAIN 

AMATRX 

CHEBY 

CMATRX 

DELTIM 

DIFF 

DJUL 

ENQ 

ERROUT 

GETTAP 

INPUT1 

MA3331 

PAGER 

READE 


SUBROUTINES AND FUNCTIONS 
NEEDED TO CREATE A SLP FILE FROM A JPL TAPE 


— identify and initiate the permanent file maintenance operations 

— computes transformation matrix which rotates from selenocentric 
to selenographic 

— fits the Chebyshev polynomial of degree mn through the points F . 
This polynomial is then converted to its equivalent power series 
and scaled to the closed interval [XMIN, XMAX]. 

— computes transformation matrix which rotates mean equator and 
equinox of 1950.0 to true-of-date 

— computes time in seconds relative to a reference date, or the 
Julian date, given a packed calender date 

— calculates the differences between any two time points in the 20 th 
Century 

— computes modified Julian date of a given Gregorian date after 
1950.0 

— (a) contains entry point — DEQ 

(b) provides data set integrity for GTDS data sets 

retrieves error messages from the permanent data base and sends 

it to the printer 

— reads and handles DE-19 tape for READE 

— controls retrieval of the ephemeris data from the JPL tape 

— computes the product of a 3 x 3 matrix and a 3 x 1 matrix 

— (a) contains entry point- PAGE NP 

(b) provides paging control to printer 

— retrieves, interpolates, and translates ephemeris data from the 
JPL Ephemeris Tape 
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RPDATO 

RYMDI 

SETDAF 

SETSLP 

SLPEPH 

SLPPF 

SLPTAP 

STADRO 

TDIF 


— obtains current date from OS/360 system 

— separates 6-digit packed calender date into two-digit words 

— supplies data control block information to FORTRAN I/O routines 

— reads options and initializes creation of the SLP file 

— controls the retrieval of the interpolated and translated JPL 
ephemeris data and computes the Chebyshev polynomial coefficients 
for output 

— (a) contains entry points: JPLTPF & SLPTPF 
(b) creates SLP ephemeris permanent file 

— generates SLP ephemeris tape 

— provides accounting information 

— computes time differences 
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APPENDIX C 


OVERALL FLOW OF THE GENERATION PROCEDURE 


This chart shows those routines in the CREATE3 program which are called when 
creating a SLP file. 


MAIN 



1 

GETTAP 
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APPENDIX D 


COMMON BLOCKS ESSENTIAL TO THE GENERATION PROCEDURE 


COMMON blocks used which transmit the data between the routines that create 
the SLP file: 

CETBL1 

CETBL3 

CETBL4 

CETBL9 

CHEV 

CONST 

FILES 

INPUT 

SAVE 

SLPOPT 

SLPREC 

SWITCH 

TIMCOF 

XLABEL 
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APPENDIX E 


ABBREVIATED GTDS SYSTEM USAGE OF THE SLP FILE 


The COMMON BLOCKS used by the routine EVAL and the routines called by 
EVAL which enable the data to be transmitted between the routines. 


SUBROUTINE 

COMMON BLOCKS 

EVAL 

CONST, FILES, FRC, 

SATPOS, SLPOPT, SLPREC, SWITCH 

ERROUT 

FILES, SWITCH 

TDIF 

TIMCOFF 

POLMOT 

TIMCOFF, SATPOS 

THETAG 

CONST, SATPOS, SLPOPT, 
SLPREC, SWITCH 
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NASA— GSFC 


£\ML Z 

COMPUTES THE 
POSITION AND 
VELOCITY VECTORS. i 
AND THE ROTATION / 
MATRICES 


ERROUT 


M 

i 

to 


RETRIEVES THE 

/appropriate error \ 

TEXT FROM THE 
PERMANENT DATA 
VBASE AND TRANSMITS/ 
\IT TO THE PRINTER j 


DELTIM 


COMPUTES THE 
TIME IN SECONDS 
RELATIVE TO A 
REFERENCE DATE. OR 
THE JULIAN DATE. 
GIVEN A PACKED 
V CALENDER OATE / 


DJUL 


COMPUTES THE 
MODIFIED JULIAN 
DATE OF A GIVEN 
GREGORIAN DATE 
AFTER 1950.0 


OBSTIM 


COMPUTES THE 
TIME ELAPSED 
FROM DAY TO 
YMD, HMS IN 
SECONDS OF 
A.1 OR UTC 




COMPUTES 

TIME 

DIFFERENCES 


THETAG 

COMPUTES THE 
GREENWICH 
HOUR ANGLE 
AT TIME 
T 


POLNOT 

COMPUTES THE 
X AND Y 
POLAR MOTION 
ANGLES AT A 
GIVEN DATE 


DJULI 


COMPUTES 
MODIFIED JULIAN 
DATE AND FRACTIONAL / 
PART OF A JULIAN 
DAY IN SECONDS 


